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Abstract 

It is shown that exact solutions may be found for the energy eigenvalue problem 
generated by the class of semirelativistic Hamiltonians of the form H = + 

V, where D is a non-local potential with a separable kernel of the form V{r,r') = 
Explicit examples in one and three dimensions are discussed, 
including the Yamaguchi and Gauss potentials. The results are used to obtain lower 
bounds for the energy of the corresponding A-boson problem, with upper bounds 
provided by the use of a Gaussian trial function. 
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1 Introduction 


We study semirelativistic problems in which the Hamiltonian H has the rel- 
ativistically correct expression K{p^) = y/m? + p = |p|, for the energy of 
a free particle of mass m and momentum p, and an added static interaction 
potential V. The Hamiltonian is therefore given by 

H = ^Jvn? + p^ + V. (1.1) 
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The eigenvalue equation Hip = Eip is usually called the spinless Salpeter 
equation [1] [2], For many potentials, this Hamiltonian can be shown [3] to be 
bounded below and essentially self-adjoint, and its spectrum can be defined 
variationally. From the point of view of solvability, these features represent 
signihcant technical advantages over the more-complete Bethe-Salpeter for¬ 
mulation. There is, however, one remaining difficulty, namely the non-locality 
of the kinetic-energy operator. 

The ‘usual’ multiplicative potential operator of elementary quantum mechan¬ 
ics is generated by a special kernel of the form V{x,x') = V{x)S{x,x'). Thus 
we have 

OO 

{Vip){x) = j V{x)6{x,x')ip{x')dx'= V{x)ip{x), (1.2) 

— OO 

and this special form makes V a local ‘multiplicative’ operator. Since (with 
h = 1) the Schrodinger kinetic-energy operator /{2m) = —S'/./{2m) is also 
local, the non-relativistic Hamiltonian is a local operator. By contrast, the 
kinetic-energy operator K = ^m^ + in the semirelativistic problem is non¬ 
local and this is the source of many of the difficulties encountered with the 
corresponding eigenvalue problem. The action of K is dehned [3] in terms of 
the Fourier transform E{ip) = ip. Thus, in one dimension, we have explicitly: 

E{kip){k) = Vm^ + k^ ip{k), (1.3) 

where 

OO OO 

'ip{k) = {'il^\k) = [ {'ilj\x)dx{x\k) = —= [ e~''^'^%lj{x)dx. (1.4) 

J 27r 

— OO —OO 

The main purpose of the present article is to show that with separable poten¬ 

tials, the non-locality of both terms in the Hamiltonian allows us to solve the 
eigenproblem exactly, up to a dehnite integral. This result, in turn, allows us 
to hnd a lower bound to the energy of the corresponding iV-boson problem in 
which the particles interact pairwise. The class of potential kernels we shall 
consider may be written (for the one-body problem in one dimension) 

n 

V{x, x') = -vY^ fi{x)gi{x'). (1.5) 

i=l 

Such potentials have been studied as models for a variety of physical problems 
[4,5,6,7,8]. Our main general results for a single particle in one and three 
dimensions are proved in Section 2. In Section 3 we look at some exponential 
examples in one dimension and in Section 4 we solve the eigenproblem for the 
three-dimensional Yamaguchi [4] and Gauss potentials. In Section 5 we apply 
the results to study a system of N identical bosons interacting pairwise in three 
dimensions via a non-local Gauss potential: the one-particle exact solutions 
provide an energy lower bound, to which we adjoin a variational upper bound 
derived with the aid of a scale-optimized Gaussian trial function. 
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2 Exact solutions 


For definiteness, we first solve the problem with one separable potential term 
in one spatial dimension. Thus we suppose that the kernel of the potential 
operator V has the form 


V(x,x') =-vf(x)g(x'), (2.1) 

where n is a postive coupling parameter. The eigen equation for the semirela- 
tivistic one-body problem becomes 

OO 

\Jm? ip{x) — J vf{x)g{x')ip{x')dx'= Eip{x). (2.2) 


If we represent the Fourier transforms by = i/j,E{f) = f, and E{g) = g, 
then Eq.(2.2) becomes 

Vm^~+^ — vcf{k) = Eip{k), (2.3) 

where the constant c is given by 


OO 

c = J g{k')'ijj{k')dk'. 

— OO 


Thus if) is given by 


If we now multiply both 
following formula relating 


^{k) = 


cvf{k) 


(2.4) 


+ W — E 

sides of (2.4) by g{k) and integrate, we hnd the 
the reciprocal coupling to the energy E\ 


1^7 f{k)g{k)dk 

V J v^m^TjTp' _ £] 

— OO 


(2.5) 


Equations (2.4) and (2.5) show that if there is a solution for given f{x) and 
g{x), then this solution is unique (up to a phase). The corresponding energy 
eigenvalue is now determined by (2.5) since n is a monotone function of E. 
It is clear that + k? > m. If we write E = E{m) = m + e{m), then, 
for bound states, e = E — m < 0. Consequently we have V — E > 0. 
That is to say, there are no mathematical singularities arising from this factor 
in the various integrands. In the large-m limit, the problem approaches the 
corresponding non-relativistic case since ^/m'^ + k"^ — m ~ k‘^/{2m); moreover, 
this approach is from below since ^/m? + k'^ — m < k"^/{2m). In the examples 
we shall consider, the function E{m) — m approaches the non-relativistic m- 
dependence as the mass increases. Another interesting special case for the 
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semirelativistic problem is the ultra-relativistic limit m —0. This, of course, 
has no natural non-relativistic counterpart. 


When there are more than one term in the separable potential we have 

n 

V{x, x') = Vifi{x)gi{x'), (2.6) 

i=l 

We now dehne the constants by 


OO 

Ci = y gi{k')ijj{k')dk', 

— OO 


(2.7) 


and the formula (2.4) for the wave function in this more-general case becomes 

n ~ 

E Vifi{k)ci 
i=l 

\/m? + k? — E 

The relation between the coupling parameters and the eigenvalue is now ex¬ 
pressed by the condition that the linear equations for the constants {cj}Ei 
are non trivial. If we dehne the matrix elements of the n x n matrix J by 


( 2 . 8 ) 


^(/c) 


Jji Xi 


OO 



— OO 


9j{k)fi{k)dk 

'Jm? + kP‘ — E' 


(2.9) 


then the more-general eigenvalue formula, corresponding to (2.5), may be writ¬ 
ten 


det(/ — J) = 0. 


( 2 . 10 ) 


There are similar results in three spatial dimensions. We consider one-term 
central potentials of the form 

V(r,r') =-n/(r)c/(r'), (2.11) 


where r = |r|. In this case the Fourier transform E{f) of /, for example, takes 
the form 

0 

in which k = |k|. Similar reasoning to that of the one-dimensional case then 
yields the solution formulae 



1 , OO 
2 \ 2 1 


TT 


sin{kr)r f {r)dr 


( 2 . 12 ) 


^{k) 


cvfjk) 

\/w? + kP' — E 


(2.13) 
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and 


= An 


f{k)g{k)k‘^dk 
^/m? + k"^ — E 


(2.14) 


These results can also easily be extended to potential kernels with a sum of 
separable terms. 


3 Problems in one dimension 


We now consider some examples. Since the general solution is given in Sec¬ 
tion 2, the purpose of the examples is to demonstrate that exact solutions are 
indeed feasible. We solve the first problem in some detail and then present 
summary solutions and results for a selection of other problems. 


3.1 The one-term exponential potential 


We consider the potential 

V{x,x') = —vf{x)f{x') = v, a > 0. 

The potential factors in momentum space are therefore given by 


(3.1) 


m = 


cos{kx) f {x)dx = 


TT Vl -I- a^/c^/ ’ 


(3.2) 


and the formulae (2.4) and (2.5) for the momentum-space wave function and 
the corresponding eigenvalue become 


^(A;) = 


cav 


{l+aV) (ym? + B-E) 


(3.4) 


and 


1 

'— 1 1 

7 Pik)dk 


1 

1 

1 

X) 

f 

1 

-\- k'^ -\- \ 

O ' 

\E\ 

n J 
0 

(1 -|- a^/c^)^ [ym? -\- k'^ — E^ 


(3.5) 


This equation may be inverted to give E for each choice of the parameter set 
{a,m,v}. In Fig.(l) we exhibit the m dependence of E — m for a = 1 and 
V = {1, 2, 3}. In the Schrodinger limit, m oo, we find 


e(m) = E(m) — m = —. 

a 


(3.6) 
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Meanwhile for the ultrarelativistic case m = 0 we have 

1 a? {2 + 2a^e^ + Saevr + a^e^vr + 4 ln(—ae)) 

V (1 + a^e^Yn 

The graphs shown in Fig.(l) are consistent with these relations. We see that 
this semirelativistic problem is indeed exactly soluble. 
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- 1.5 

E - m 

-2 

- 2.5 

-3 

- 3.5 

0 1 2 3 4 5 

m 

Fig. 1. Plots of exact semirelativistic energies E{m) — m for the nonlocal exponential 
potential V{x,x') = —I for three values of the coupling v. 


3.2 A two-term exponential potential 

We consider now the case 

V(a;, x') = — Vb, a, b >0. (3.8) 

In particular, if we choose the explicit values a = 1, 6 = 2, = u;, = 1, the 

secular equation (2.10) becomes 

(1 — Jii)(l — J22) — J12 = 0, (3-9) 
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where the integrals are given by 


4 7 dk 

71 + - E)' 


(3.10a) 


and 


16 7 dk 

TT 7 (1 + Ak'^Y{7+ k"^ — Ey 


OO 

Sr dk 

(1 + A;2)(l + 4A;2) {y/m^ + k^ - E)' 


(3.106) 


(3.10c) 


Thus for m = {0, 0.5, 1} we hnd respectively from (3.9) that E = {—1.14462, — 
0.814543, —0.36131}. In Fig.(2) we exhibit a graph showing E — m as a func¬ 
tion of m for this problem. 



Fig. 2. Plots of exact semirelativistic energies E{m) — m for the nonlocal two-term 
exponential potential V{x,x') = —I). 
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4 Problems in three dimensions 


The Yamaguchi potential [4] has a potential kernel given by 


V(r, r') = —V 


e-Pr 


r 



V, /9 > 0. 


(4.1) 


Because of the volume measure r‘^dr in three dimensions cancels the singular¬ 
ities in the Yukawa-type factors, this problem is very similar to the exponen¬ 
tial potential in one dimension. The wave function and eigenvalue formula are 
found from (2.13) and (2.14) to be respectively 


^{k) = 


(A;2-f/32) + k^-E) 


(4.2) 


and 


1 

V 


k‘^dk 


(k^ + iPf (VSPTF-B 


(4,3) 


Thus the energy E may be found from (4.3) as a function of the positive 
parameters {m,/5,n}. 


Similarly, for the Gauss potential we have the kernel 

V(r,r') = —^ n, /3 > 0. (4.4) 


The corresponding wave function and eigenvalue formula in this case are given 


by 

ipik) = =- 

^ ^ Vm^ + k^-E 

(4.5) 

and 

1 Att J e~^^^^k'^dk 

(4,6) 


V f3^ J y/m? + k‘^-E 


5 The semirelativistic Y-boson problem 


In this section we consider a system of N identical bosons interacting pairwise 
in three spatial dimensions. The Hamiltonian for the system may be written 

N I N 

H = Y.(m" + Pp^+ E Y 

2=1 Jf>2 = l 

8 


(5,1) 



where, for a single particle, the action of the Gauss potential is given by 

G'^(r) = —V J , v > 0. (5.2) 

We shall consider two distinct approaches. First we obtain a lower bound to 
the lowest iV-body energy E with the aid of a scaled one-body problem and 
secondly we hnd an upper bound with the aid of an iV-body Gaussian trial 
wave function. 


5.1 The lower bound 


If we suppose that 'h is the exact (unknown) iV-boson wave function, then 
boson symmetry implies that E = (d', H'$) = (\h, /i\h), where h is a two-body 
Hamiltonian given by 


h 


N 

Y 


[m^ + pj)2 + + P2)2 + (iV - l)Ci2 


(5.3). 


If new coordinates for the two-body problem are r = ri — r 2 and p = ri -|- r 2 , 
then the individual momenta are given by tt ± p, and, by using the lemma 
of Ref. [11] to ‘remove’ the operator vr from within expectation values, the 
two-body operator h may be replaced by TC, where 


n 


N 


{m^ + p^)2 + UN - l)V 


Thus we conclude that 


(5.4) 


B = (T,//I') = (>I',W4') > f = Bi, (5.5) 


where S is the bottom of the spectrum of the one-body operator Ti. By com¬ 
paring (5.4) with (4.6) we see that 

1 2tt °r 

{N-l)v ^ ^ y - EJN' ^ 

Thus, for each choice of the parameters m and (3, (5.5a) implies that E^/N is 
a function of v{N — 1). In the special case m = (3 = 1, we write this function 
as fi so that we have 

EL/N = fLiviN-l)). (5.56) 

We note the special critical coupling Uc dehned by f{uc) = 0 is given by 
y(iV- l)|c = 0.527485. 
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0 1 2 3 4 5 


(N-1) V 


Fig. 3. Upper and lower bounds for E/N against {N — l)v for N identical bosons 
interacting via the Gauss pair potential V{r,r') = ^ . The lower curve 

(L) is for N = 2 and the upper curve (U) is for N = oo; the curves for 2 < N < oo 
lie between these two. 

5.2 The upper bound 


For a variational upper bound we adopt explicit relative coordinates. Jacobi 
coordinates may be defined with the aid of an orthogonal matrix B relating 
the column vectors of the new [pi] and old [r*] coordinates given by [p*] = B[ri]. 
The first row of B defines a center-of-mass variable pi with every entry 1/v^, 
the second row defines a pair distance p 2 = (ri — r 2 )/v^, and the Ath row, 
k > 2, has the hrst k — 1 entries B^i = 1/^k{k — 1), the Ath entry Bkk = 

~\l ~ 1)/^) ^he remaining entries zero. We dehne the corresponding 
momentum variables by [tt*] = (i?“^)*[pj] = B\pj\. The trial wave function we 
use is given by 


<F(P 2 , P3, • • •, Pn) = Ce 2 S -2 = G n 0(Pi) 

i =2 


(5.6) 
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where a > 0 and C* is a normalization constant. This function is symmetric 
in the individual position coordinates and also in the iV — 1 relative 

coordinates meanwhile it has the unique factoring property shown. 

These facts enable us [11] to express the expectation of the full Hamiltonian 
H in the form 


E < (T, ifT) = {{m^ + 2Ap2)^ + {N - 1)H)0), (5.7) 


where A = (iV —l)/iV, the potential operator V has the Gauss kernel (4.4), and 
a is to be used as a variational parameter. We therefore obtain the following 
expression for the upper bound Eu in the special case m = {3 = 1 


N - N 



min 

s>0 


9{s^) 


— SUTT^ 


( 2 As ^)2 
(l + 4As2)3 ’ 


(5.8a) 


where the monotone function g is given by 


9{x) 



2x + t‘^ 


t^dt 


xe^Ki{x). 


In this last expression, Ky{x) is a modihed Bessel function of the second kind 
[12]. The result of the minimization in (5.8a) yields Ejj/N as a function fu of 
(TV — l)n and A, where | < A < 1. We have 


EulN = fu{v{N-l),\). 


(5.85) 


Thus we obtain a different upper-bound curve for each A = (TV — 1)/TV. These 
curves do not intersect. In Fig.(3) we exhibit the lower curve fiiv^N — 1)), 
valid for all iV, the upper curve fuiv{N — 1), |), for N = 2, and the upper 
curve fu{v{N — 1),1), for N = cx). For the case N = 2 the general lower 
(all N > 2) and particular upper bounds (TV = 2) are so close that they are 
indistinguishable on the graph; we have, for example, /l( 1) = —2.56844 and 
/c/(l, |) = —2.5651 approximately. Thus the scale-optimized Gaussian trial 
function is very effective for all N, and particularly so for N = 2. The apparent 
straightness of the energy curves can perhaps be understood by reasoning such 
as the following: for the lower bound (5.5a), the Gaussian in the integrand 
decays rapidly to zero, thus the mean-value theorem tells ns, for a given n, 
that {N — l)n = A — B{El/N)-, it remains, of course, to explain why A and B 
vary very slowly with n.However, with exact analytical results available (for 
both bounds), we do not have to look for more analytical approximations. 
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6 Conclusion 


We have shown that exact solutions can be found to semirelativistic eigenvalue 
problems when the potential has a kernel that is a sum of separable terms. This 
immediately extends, of course, to the wider class of kernels. It may be pos¬ 
sible to use such exact solutions to approximate the spectra generated by local 
potentials. The non-relativistic many-body problem with non-local potentials 
has already been studied [9] and the present paper extends these results to the 
corresponding semirelativistic case. We have obtained tight bounds for the lo¬ 
cal semirelativistic iV-body problem with local harmonic-oscillator potentials 
V{r) = nr^, and somewhat weaker bounds for convex transformations g{r^) 
of the oscillator [10]. The work reported in the present paper will no doubt 
help us to extend these semirelativistic many-body results to wider classes 
of potentials. It is very helpful when the lower bound itself, which is derived 
from a scaled one-body problem, can be found exactly. Improvements in the 
general lower bound await a treatment based on Jacobi relative coordinates; 
this has already been achieved in particular for the oscillator; the search for 
an improved general lower bound can now beneht from a non-oscillator test 
model for which there is also an accurate variational upper bound. 
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